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We investigate the properties of a single photon generated by a solid-state emitter sub- 
ject to strong pure dephasing. We employ a model in which all the elements of the 
system, including the propagating fields, are treated quantum mechanically. We analyt- 
ically derive the density matrix of the emitted photon, which contains full information 
about the photon, such as its pulse profile, power spectrum, and purity. We visualize 
these analytical results using realistic parameters and reveal the conditions for maxi- 
mizing the purity of generated photons. 
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1. Introduction 

Cavity quantum electrodynamics (QED) is currently one of the hottest research 
topics in atomic physics. In particular, quantum-mechanical interactions between atoms 
and cavity photons have been intensively studied. 1 Solid-state cavity QED systems 
composed of semiconductor quantum dots and cavities have recently been attracting 
much attention since they are suitable for creating compact optical devices. 2 Strong 
coupling between a single dot and a cavity has been confirmed through a large vacuum 
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Rabi splitting. 3-5 This strong coupling has been applied to fabricate a single-dot laser 6 
and to generate nonclassical light including single photons. 7-11 Excellent performances 
have been reported in generating indistinguishable photons 12 ' 13 and entangled photon 
pairs, 14 both of which are useful for quantum information processing. 15 

In contrast with real atoms, semiconductor quantum dots are strongly influenced by 
environmental noise sources such as phonons and background carriers. The fluorescence 
spectra of solid-state and atomic cavity QED are qualitatively different. When a solid- 
state system is excited by a pump light of the dot frequency, a spectral peak appears at 
the cavity frequency in spite of the large detuning between them. 3-5 This phenomenon 
has not been reported in atomic cavity QED systems. Further experimental studies have 
characterized this peak more thoroughly 16-22 and have revealed that the fluorescence 
at the cavity frequency is due to radiative decay of the dot. Subsequent theoretical 
studies accounted for pure dephasing of the dot through the stochastic Schrodinger 
equation 23,24 or the Master equation 25-27 and successfully explained the peak at the 
cavity frequency. The influence of pure dephasing on radiative decay of the dot can be 
understood in terms of the quantum Zeno and anti-Zeno effects. 25 ' 28 

Therefore, when designing a single-photon source using solid-state cavity QED sys- 
tems, it is crucial to quantitatively consider pure dephasing of the dot. The performance 
of such photon sources should be evaluated from two aspects. One is the collection effi- 
ciency; namely, the probability that the emitted photon is transferred to the intended 
spatial mode (i.e., the radiation pattern of the cavity). This has been discussed in sev- 
eral studies in terms of the ratio of radiative decay rates. 23 ' 26 ' 27 The other one is the 
indistinguishability of generated photons, which can be measured by two-photon in- 
terference experiments and is evaluated by the purity. Single photons with high purity 
are required for quantum information processing, particularly for constructing scalable 
quantum circuits. 

In this study, we investigate the properties of a single photon emitted by a solid- 
state cavity QED system and quantitatively observe the effects of pure dephasing. We 
treat the five elements of the overall system (the dot, the cavity, radiation leaking from 
the cavity, non-cavity radiation modes, and the environment causing pure dephasing 
of the dot) as active quantum-mechanical degrees of freedom, and analytically derive 
the density matrix of the emitted photon in the real-space representation. This density 
matrix contains full information about the emitted photon, including its pulse profile, 
the frequency spectrum, and the purity. These quantities are observed as functions of 
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Fig. 1. Schematic illustration of the solid-state cavity QED system considered. It consists of a quan- 
tum dot, a cavity, photon leakage from the cavity (b field), non-cavity radiation modes (c field), and a 
reservoir field, which causes pure dephasing of the dot (d field). 



the pure dephasing rate of the dot. We reveal the optimum condition for maximizing 
the purity of the emitted photon. 

2. Model 

We investigate radiative decay of an excited quantum dot placed inside a cavity, 
as illustrated in Fig. 1. This solid-state cavity QED system consists of the following 
five components: (i) a quantum dot, (ii) a cavity, (iii) a photon field leaking from the 
cavity (referred to as the b field hereafter), (iv) non-cavity radiation modes (c field), and 
(v) a reservoir field, which causes pure dephasing of the dot (d field). 34 The annihilation 
operators corresponding to these components are respectively denoted by a, a, b k , c k , 
and d k , where k is a one-dimensional wavenumber. Note that a is a Pauli operator, 
whereas the other operators are bosonic. Putting h = c = 1, the Hamiltonian of the 
overall system is given by 

n = n + n 1 + 'H2 + n 3 , (i) 

Ho = Ud.o~^o~ + uj c a)a + g(a^a + aV), (2) 



Hi 



n 2 



dk 



dk 



dk 



kb\b k + V^/(27r)(a f 6 fc + b\a) 
kc\c k + v / 7/(2vr)(o- t c fc + c\a) 
kd\d k + y / 7~/W tcr ( c 4 + d k 



(3) 



(4) 



(5) 



J. Phys. Soc. Jpn. 



DRAFT 



The parameters are defined as follows (see Fig. 1). cu d and co c respectively denote the 
resonance frequencies of the dot and the cavity, g represents the coupling between 
them, k is the escape rate of cavity photons, 7 is the radiative decay rate of the dot 
into non-cavity modes, and 7 P is the pure dephasing rate of the dot. describes the 
Rabi oscillation between the dot and the cavity (Jaynes-Cummings Hamiltonian), "Hi 
describes leakage of a cavity photon to its radiation pattern, "H 2 describes the radiative 
decay of the dot into unintended directions, and "H 3 describes pure dephasing of the 
dot. We can confirm that Af = <tV + a)a + f drb\b r + / drclc r commutes with the 
Hamiltonian. Therefore, the number of excitations is conserved in the dot, cavity, b and 
c fields. 

We assume that the dot is initially (t = 0) in the excited state while the other fields 
are in their vacuum states. Then, denoting the overall vacuum state by |0), the initial 
state vector is given by 



The Hamiltonian of eq. (1) and the initial state vector of eq. (6) form the basis of our 
analysis. 

For later convenience, we introduce the real-space representation of the b field (cavity 
leakage). It is defined by 



The r < (r > 0) region represents the incoming (outgoing) field. c r and d r . can be 
formally defined in a similar manner. Our main concern lies in the properties of single 
photon emitted in the b field. 

To model the pure dephasing of the dot, the d field interacts with the dot so as 
to conserve the dot excitation. Using eqs. (2) and (5), the dot Hamiltonian can be 
rewritten as [ui d + f(t)]a^a, where f(t) = y/2^ p [d (t) + dl(t)} is the fluctuation of the 
dot resonance frequency induced by the d field. Using eqs. (6) and (10), we can confirm 
that (f(t)f(t'))i = 27 p (5(t — t'), where (• • -)j = ■ ■ ■ Therefore, the present model 
assumes a white noise spectrum for fluctuation of the dot resonance. 

3. Analysis 

In this section, we present analytical results that are rigorously derivable from the 
model of §2. We solve the time evolution of the overall system within the input-output 
formalism 29 ' 30 and derive several formulae to characterize the emitted single photon 



m = ^|o>. 



(6) 




(7) 
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(density matrix, pulse shape, spectrum, and purity). These analytical results are visu- 
alized in the next section for specific parameters. 

3.1 Heisenberg equations 

Here, we present the Heisenberg equations for the system (er, a) and field (&&, dk) 
operators. Deriving the raw Heisenberg equations for the field operators from eq. (1) 
and transforming them into real-space representations, we obtain the following relations 
that connect the incoming (r < 0) and outgoing (r > 0) fields: 

b r (t) = 6 r _ t (0) - iy/K9{r)9(t - r)a(t - r), (8) 

c r (t) = cr- t (0) - i^e{r)6{t - r)a(t - r), (9) 

d r (t) = d r _ t (0) - i^/¥f p Q(r)e(t - r)a ] (t - r)a(t - r), (10) 

where 6(x) is the Heaviside step function. From the raw Heisenberg equations for the 
system operators and the above input-output relations, the Heisenberg equations for 
system operators are given by 

-^-<7 = —iuidcr — ig(l — 2o~^o~)a — i(l — 2a^a)N c (t) — % N ] d (t)a + a*N d (t)\ , (11) 



d_ 
It 



—iuj c a — igcr — iN b (t), 



(12) 



where Ud = ojd — i{j /2 + ^ p ) and u c = u c — in/ 2 are respectively the complex frequencies 
of the dot and the cavity, and the noise operators are defined by Nj,(t) = y/K,b-t(0), 
N c (t) = ^70-4(0), and Nd(t) = y/2j p d- t (0)- Note that the noise operators are the 
initial-time operators and consequently Nj(t)\0) — (j — b,c,d). 



3.2 State vector 

The state vector of the overall system at an arbitrary time t is determined by 
\ip(t)) = e~ Mt \ip(0)). Since the initial dot excitation is conserved in the dot, cavity, 6, 
and c fields, the state vector can be written as 

a (t)a ] + (3 Q {t)a ] + / rfr 7o (r,t)6j+ / dr5 (r,t)cl 



W)) 



m=l 



X 



t)a ] + (3 m (x,t)a i + I dr~f m (r,x,t)bl+ I dr5 m (r,x,t)$ 



<---4 m |o), 



(13) 

where m denotes the number of excitations in the d field and J d m x denotes a multi- 
dimensional integral with respect to x = (x±,X2, ■ ■ ■ ,x m ). We can set x 1 < ■ ■ ■ < x m 
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Fig. 2. Ai and A2 in the complex plane. Ai 
their traces as g increases. 



-iuid and A2 = — iuj c when g = 0. Dotted lines show 



(14) 



without loss of generality. As we show later, these coefficients are nonzero only when 
< r < xi < • • • < x m < t. 

First, we discuss ao and From eq. (13), we can confirm that ao(t) = (o"(t)o^(0)) 
and (3 (t) = (cOty^O)), where (• • •) = (0| • • • |0). From eqs. (11) and (12), their equa- 
tions of motion are given by 

d_ I <*o(t)\ _ l-itid -ig \ I a (t) 
Jt \Po(t) )~\-ig -&c) \Po(t) / 
with the initial conditions a (0) = 1 and (3q(0) = 0. The solutions are given by 

a (t) = Aie Xlt + A 2 e X2t , (15) 

fa{t) = B 1 e x « + B 2 e x *, (16) 

where Ai and A2 are the two eigenvalues of the 2x2 matrix in eq. (14) (see Fig. 2), 
A 1 = (Ai+iC c )/(Ai-A 2 ), A 2 = (A 2 +i£5 c )/(A 2 -Ai), and B 1 = -B 2 = -ig/(\i-\ 2 ). The 
real parts of Ai and \ 2 are always negative and consequently ao and (3o vanish as t — > 00. 
Equation (13) also implies that 70 (r, t) = (b r (t)a^(0)} and 5 (r,t) = (c r (t)a l ~(0)). From 
eqs. (8) and (9), we have 

-Yo(r,t) = -iy/RPo(t-r), (17) 

5 (r,t) = -iyfya (t-r). (18) 

Rigorously, 9(r)6(t — r) should appear on the right-hand sides of these equations. How- 
ever, below, we implicitly assume 0<r<a;i<---<a; m <t and omit the Heaviside 
functions. 
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Next, we proceed to investigate higher-order quantities including reservoir (d field) 
excitations. We consider a.\ and j3i as examples. Applying the same reasoning as that 
for 7q and 8 , we have ai(x,t) = — i^2 r y p {a(t)a^ (t — x)a(t — x)crt(O)) and f3 1 (x,t) = 
—i ■ s j2^ p {a{t)a\t — x)o~(t — x)o~'(0)). Thus, we need to evaluate the three-time corre- 
lation functions, ((T(t 2 )(T* (t 1) a (h) (7^(0)) and (a(t 2 )(T*(ti)(T(ti)(ji(0)) with t 2 > h > 0. 
Their equations of motion with respect to t 2 have the same form as eq. (14) and their 
initial values {t 2 —> ti) are a (ti) and 0, respectively. This implies that the two-time 
correlation functions can be factorized as (o"(t 2 )o^(£iV(£i)<^(0)) = cto(t 2 — ti)«o(^i) 
and (a(t 2 )a^ (ti)a{ti)a^ '(0)) = ftv(t 2 — ti)ct (ti). Thus, we have 

ai(x,t) = —iy^2 / y p a (t — x)a (x), (19) 

Pi(x, t) = -t^2T P a (t - x)0o(x). (20) 

Repeating the same arguments, all coefficients can be written as products of a and (3 , 
as follows: 

a m (x,t) = (-i^/2r^) m K m (x,t)a (x 1 ), (21) 
f3 m (x, t) = (-iy^~ p ) m K m (x, f)A,(:ri), (22) 
7 m (r, x, t) = (-iy/n) (-i \[2n P ) m K m (x, t)(5 Q {x l - r), (23) 
5 m (r, x, t) = (-iy/i) (-iy^) m K m (x, t)a {x 1 - r), (24) 



where 



K m (x, t) = a (t - x m )ao(x m - z m _i) • • • a {x 2 - xi). (25) 



3.3 Decay of dot excitation 

The state vector of eq. (13) fully describes the dynamics of the overall system, 
including both its transient and asymptotic behaviors. In this section, as an example 
of a transient phenomenon, we analyze the decay of the dot excitation. The survival 
probability of the dot excitation is defined as P(t) = {ip(t)\a^a\ip(t)). From eqs. (13) 
and (21), we have 

P(t) = \a^{t)\ 2 + 2 lp J dx\a (t - x)\ 2 \a (x)\ 2 + ■ ■ ■ . (26) 

We here introduce the Laplace transform of |a | 2 , which is defined by C\ ao \2(z) = 
Jq 00 dte~ zt \a (t)\ 2 . It is given by 

4oop(*)= e - A : K x^ (27) 

m,n=l,2 " 
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Fig. 3. /Zj (j = 1, • • • ,4) in the complex plane. When j p is absent, \x\ = Ai + AJ, /12 = Ai + Aj, 
= ^2 + AJ, and ^4 = A 2 + A 2 . Dotted lines indicate their traces as 7 P increases. Real parts of /ij are 
always negative for any -f p . 



where A ii2 and A 1:2 are denned in §3.2. The Laplace transform of P(t) is then given by 

Cp(z) = C ":f Z \ z y (28) 

P(t) is obtained by analyzing the poles of this function in the z-plane. We denote the 
four roots of the equation 1 — 2j p £\ ao \2(z) = by fij (j — 1, • • • ,4) (see Fig. 3). P(t) is 
then given by 

4 

P(t) = J2E^ t , (29) 
i=i 

E = Y\ m ',n' =1,2^0 ~ X rn> ~ A*,) A m A* n 

Note that the real parts of jij are always negative and that the survival probability P(t) 
vanishes in the t — > 00 limit, as expected. 

3.4 Density matrix of emitted photon 

In the t — > 00 limit, the initial dot excitation is completely transformed into a 
photon propagating in the intended mode (b field) or in unintended directions (c field). 
In the following subsections, we analyze the photon emitted in the 6-field. It is fully 
characterized by its density matrix p(t). In the real space representation, the matrix 
element p(r, r', t) is given by 

p(ry,t) = (m\b%\m)- (31) 
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We make the following three comments regarding this quantity: (i) p(t) is Hermitian, 
namely p(r f ,r,t) = p*(r,r',t). Therefore, we need consider only the r < r' region in 
this subsection, (ii) As we will see later, p(r, r', t) = p(r — t,r' — t) in the t — > oo limit. 
This reflects translational motion of the emitted photon, (iii) Trp(t) = J drp(r, r, t) 
represents the probability of finding the emitted photon in the b field. This is unity 
when 7 = 0. 

Using eqs. (13) and (23), the matrix element can be rewritten as 
p(r,r',t) = K^{t-r)^{t-r') + 2 lp K J dx^ (x - r)/3*(rr - r')\a (t - x)\ 2 + • • • . 

(32) 

We here introduce the Laplace transform of /3o/3q, which is defined by Cp p*(r, r', z) = 
J °° dte~ zt /3 (t - r)/3*(t - r'). It is given by 

C m (r,r',z)= £ B ™ K e^'-^'z, (33) 

m,n=l,2 Z - X rn~K 

where Ai^ and B\2 are defined in §3.2. The Laplace transform of p(r, r', t) is then given 
by 



«' y '"=r^Ti' (34) 

1 - 2^ p L\ ao \2{z) 



By analyzing the poles of this function in the z-plane, p(r, r', t) is obtained as follows: 
p(r, r>, t) = J2 E ft^^^), (35) 

j=l m=l 

ri m ',n'=i,2(/ i i ~~ Am' - A*,) B m B* 

pjm = fr — u- — 1T\ ^ u - a — a* ' (36) 

where < r < r' < t and pj (j — 1, • • • , 4) are defined in §3.3. We can check that this 
quantity depends on only two variables, r — t and r' — t. 

3.5 Pulse shape 

The pulse shape of the emitted photon is characterized by the intensity distribution 
/(r, t) = {ip(t)\blb r \i/j(t)) = p(r, r, i), which is the diagonal element of the density matrix. 
This quantity is real and positive in < r < t. By setting r' = r in eq. (35), we have 

4 

/(r,0=E^ eW( '" r) ' ( 3? ) 

3=1 

, Hm',n'=l,2(/ i j ~ Am' ~ A*,) -B m -B* 

" ' rW^"-*) A, ; , A;, 
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3. 6 Frequency Spectrum 

The frequency spectrum of the emitted photon is defined by S(k, t) = 
(i J (^)rk^k\'4 ) ('t))- Apparently, S(k,t) becomes independent of t in the t — > oo limit, 
and we are interested in S(k) = \im t ^ ^ S(k,t). By definition, S(k,t) is the Fourier 
transform of the density matrix element: 

S(k,t) = ^- J J drdr'e ik{r '- r) p(r,r',t). (39) 

We consider the Laplace transform of S(k,t) defined by Cs(k,z) = / °° dte~ zt S(k,t). 
Using eqs. (34) and (39), it is given by 

£ S (M) = ~ C ' m ?' Z) ,„ (40) 
SV ' 27rl-2 7p £| Q0 | 2 (z)' 1 ; 



where 



dtJJ drdr'e ik ^-^- zt p*(t-r')(3*(t- 



oo 

-zt 



/ dte 
Jo 



2 

£-%(e---e^) 



(41) 



Note that ^^ o/3 » has a pole at z — 0. Since our interest lies in the t — > oo limit of 
S(k,t), we need to investigate the pole of Cs{k,z) at z = only. Therefore, S(k) = 
(k/27t)[1 — 27 P £| Q , p(0)] _1 Res^=o[£|a 0/ g*(A;, z)]. After some calculations, we obtain 

S ^ = Wl 212' ( 42 ) 

\(k - u d )(k - u c ) - g l Y 
where M = (ng 2 /27r)[l — 27 P £| QO |2(0)] _1 is a factor that is independent of k. This 
spectral shape was predicted by Glauber, 33 and it was confirmed in recent theoretical 
studies based on the quantum Langevin equations 23,24 and the Master equations. 25-27 
The relation to the Master equation approach is summarized in Appendix A. 

3.7 Purity 

Quantum information processing requires high indistinguishability between single 
photons. A popular measure of indistinguishability of photons is the coincidence proba- 
bility P co in two-photon interference experiments (see Fig. 4). Two solid-state emitters 
simultaneously emit single photons into two input ports of a beamsplitter. These two 
photons are mixed by the beamsplitter and are counted by photo detectors.. When 
two indistinguishable photons are simultaneously input to a beamsplitter, they always 
appear at the same output port (Hong-Ou-Mandel interference), namely, P co = 0. 
However, pure dephasing generates entanglement between an emitted photon and the 
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Fig. 4. Schematic illustration of two-photon interference experiment. The coincidence probability 
vanishes when the two input photons are completely indistinguishable. 



environment of its source, and it makes two-photon interference imperfect (P co > 0). 
The coincidence probability is related to the purity V of a photon by P co = (1 —V)/2 
(see Appendix C for the derivation). 

The purity is defined in terms of the density matrix p(t) by V = Tr[p 2 (t)]. As 
expected, this quantity becomes independent of t when t is sufficiently large. Using the 
real-space matrix element, the purity can be rewritten as 

V = lim / drdr' p(r,r' ,t)p(r' ,r,t) = lim / drdr'\p(r, r', t)\ 2 . (43) 

t— voo J t— >oo J 

Using eq. (36), we have 

3,j —1 m,m =1 j 

4. Numerical results 

The analytical results derived in the previous section are rigorous and applicable to 
any set of parameters, (to a — u c , g, k, 7, 7 P ). In this section, we visualize these results by 
employing specific parameters. Throughout this section, we assume for simplicity that 
radiative decay to unintended modes does not occur (7 = 0). 

4-1 Zeno and Anti-Zeno effects 

First, we observe the effects of pure dephasing on decay of the dot excitation. We fo- 
cus on the weak-coupling regime (k = 6g) in this subsection, where the dot decays mono- 
tonically without revival and obeys the exponential decay law to a high accuracy. The 
decay rate of the dot is well defined in this case and is given by T = lim^oo [— log P(t) / 1] . 
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Fig. 5. Survival probability S(t) for (a) resonant (uid — oj c ) and (b) detuned (ujd — u c = 600 fJ.eV) 
cases, g = 25 /xeV and n = 150 /xeV. The values of 7 P are indicated in the figures. 
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Fig. 6. Dependence of the radiative decay rate f on the pure-dephasing rate 7 p . The radiative decay 
rate is normalized by the free decay rate (i.e., the rate for 7 P = 0). g = 25 /zeV and k = 150 /ieV. 
— lj c = (solid), 66 /xeV (dotted) and 150 /ieV (dashed). 



This reduces to mirij |Re fij\, where fij is defined in §3.3. Figure 5 shows the temporal 
behavior of the survival probability P(t). In Fig. 5(a), where the dot is in resonance 
with the cavity (u>d = u> c ), the decay becomes slower as pure dephasing increases. In 
contrast, in Fig. 5(b), where the dot is detuned from the cavity {uj d — u) c = 600 /xeV), the 
decay becomes faster under small pure dephasing (7 P = 200 /xeV), whereas the decay 
becomes slower under larger pure dephasing (7^ = 12800 fieV). 

Measurements to check the survival of the dot excitation destroy the quantum co- 
herence between the excited and ground states preserving the population of these two 
states. Therefore, pure dephasing makes the same effects on the dot as continuous mea- 
surements, if measurement results are unquestioned. The changes of the decay rates 
induced by pure dephasing are often interpreted to as the quantum Zeno and anti-Zeno 



J. Phys. Soc. Jpn. 



DRAFT 



(a) 


= 50[^VJ^^^ 






7p = 200 [/ieV] 





-300 -250 -200 -150 -100 -50 

r - t [ps] 





0.03 


> 






0.02 






to 






0.01 




°1 



(b) 


I ,. o 

; \ 7 P - 30 /icV 
A-V 200|/ioV 











-50 50 

k - w c \jjeV] 




125 250 

7 P b»V] 
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effects. ' ' Previous analyses predict that the anti-Zeno effect can be observed only 
when the dot-cavity detuning is large and pure dephasing is small. 28 This agrees with 
the present numerical results. Figure 6 shows the dependence of the radiative decay 
rate T on the pure- dephasing rate 7 p . To clearly observe the Zeno and anti-Zeno effects, 
r(7 P ) is normalized by the free decay rate of T(0): r(7 P )/r(0) < 1 indicates the Zeno 
effect, whereas r(7 P )/r(0) > 1 exhibits the anti-Zeno effect. Non-monotonous behavior 
of r(7 P ) is clearly observed for large detuning. 

4-2 Pulse shape and spectrum 

In this subsection, we examine the pulse shape and the frequency spectrum of the 
emitted photon using the same parameters as those used in the previous subsection. 
First, we observe the results for the resonant (ud = uj c ) case. The pulse shapes /(r, t) of 
the emitted photon are shown in Fig. 7(a) for three pure dephasing rates j p . Each pulse 
shape is normalized [J^_ drf(r, t) = 1] since 7 = is assumed here and a single photon 
is necessarily generated in the b field. The pulse becomes longer as pure dephasing 
increases. This is consistent with the quantum Zeno effect discussed in the previous 
subsection: In the resonant case, the decay of the dot becomes monotonously slower 
with increasing pure dephasing. Figure 7(b) shows the spectra S(k) of the photon for 
the same parameters as Fig. 7(a). These spectra have a single peak at k — Ud(= oo c ) 
and are normalized [J™ dkS(k) = 1]. The spectrum broadens with increasing pure 
dephasing. This is confirmed by Fig. 7(c) in which the spectral width (defined as A = 
[J dk{k — ujd) 2 S{k)] 1 ^ 2 ) is plotted as a function of the pure dephasing rate. Thus, the 
pulse broadens in both real and frequency spaces with increasing 7 P and it is thus not 
Fourier-limited. This implies that the emitted photon is in a mixed state when 7 P ^ 0. 
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Fig. 8. (a) Pulse shape f(r,i) and (b) spectrum S(k) of emitted photon, g = 25 £teV, fc = 150 /xeV 
and ujd — uj c — 600 /t/eV. 7 P = (solid), 50 [ieV (dotted), and 200 fieV (dashed). 



The purity of the photon is discussed later in §4.3. 

Next, we observe the results for the detuned case. Figure 8 (a) shows the pulse shapes 
of the photon. The pulse shape is approximately exponential except for the oscillatory 
behavior at the very initial stage. The pulse length is inversely proportional to the 
decay rate of the dot. Figure 8(b) shows the photon spectra for the same parameters as 
Fig. 8(a). A notable difference from the resonant cases is that the spectra are doubly 
peaked with peaks at both the dot frequency Ud and the cavity frequency u c . The widths 
of these peaks are determined by [IiruUdj — 7/2 + 7 P and |Imc5 c | = tz/2. Therefore, the 
width of the peak at uj d is sensitive to pure dephasing. When pure dephasing is weak 
(7p *C k) as in atomic cavity QED systems, the dominant peak of S(k) appears at 
u d . In contrast, when pure dephasing is strong (j p ^> k) as in solid-state systems, the 
dominant peak of S(k) appears at u c . This behavior can be seen more clearly in a 
approximate form of the spectrum 

S(k) = -^SU^ + T^—S^k), (45) 

bcav[k) = (Jfc- Wc ) 2 + («/2) 2 ' { } 

S dot (k) = 77 2 , (47) 

(k - u d ) 2 + 7^ 

which is valid for \ud — oj c \ ^> g. From this expression, it is proved that the cavity 
spectrum S C3jV is dominant for 7 P ^> n, whereas the dot spectrum Sdot is dominant for 
7 P <C K. This result partly explains the detuned peaks observed in the resonance fluores- 
cence spectrum in solid-state cavity QED systems as discussed in previous theoretical 
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works. 23-25 

The mean energy of the emitted photon is evaluated by E p = J dk kS(k). Due to 
energy conservation, the mean photon energy is expected to always be identical to the 
dot frequency (i.e., E p = Ud). However, Fig. 8(b) clearly shows that the mean photon 
energy is sensitive to pure dephasing and may deviate from the dot frequency when 
pure dephasing is present. This discrepancy can be resolved by considering the energy 
released to the environment during decay, which is given by E e = J dk k(d\d k ). It can 
be shown that (see Appendix D for derivation) 

E P = ~^-^c + (48) 
k + 27 P k + 2"f p 

E e = ^—(uj d -uj c ). (49) 

K + 2^ p 

Thus, energy conservation is satisfied when the environmental energy is included (E p + 
E e = ojd)- It should be noted that while pure dephasing coupling never induces a 
transition in the dot, this coupling enables energy exchange between the dot and the 
environment. When the cavity frequency exceeds the dot frequency, eq. (49) indicates 
that the dot may absorb environmental energy during decay. 34 



4.3 Purity 

As discussed in §3.7, the coincidence probability in two-photon interference P co , 
which is a popular measure of indistinguishability, is related to the purity as P co = 
(1 — V)/2. In this subsection, we show numerical results on the purity V for various pa- 
rameter regions. Qualitative features in some limiting cases are discussed in Appendix B 
by perturbation calculation. 

We first consider the resonant case (uu c = uj). In Fig. 9 (a), the purity V is plotted 
as a function of k and j p . The purity becomes large for 7 P < min(ft/2, 2g 2 /n) (the region 
below 'Line 1' in the figure). From the figure, we can see that for generation of coherent 
photons, pure dephasing should be reduced as much as possible, and that for the fixed 
pure dephasing, k should be taken near the critical damping (k = 2g). In addition to 
this main region appropriate for highly-coherent photon generation, there is a line on 
which the purity has a small peak (denoted with 'Line 2' in the figure) in the region of 
7p ^ 9- From perturbation calculation (see Appendix B), we can show that the purity 
takes a maximum value V nmx = 3 — 2\[2 ~ 0.17 at 7 P = 2\^2g 2 /n. 

Next, we consider the largely detuned case (\ud — ou c \ ^> g). In Fig. 9 (b), the purity 
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Fig. 9. The purity calculated as a function of k/ g and 7 P /g for (a) the resonant case {uid = w c ) and 
(b) the detuned case (w f j — uj c = 8g). 



is shown as a function of k and 7 P for — u c = 8g. The purity has a large value 
below the line 7 P < min(eK/2,2g 2 / k) (the region below 'Line 1' in the figure), where 
e = g 2 /(uid — u c ) 2 is a small parameter dependent on the detuning. In Fig. 9 (b), we 
can see that the purity for the detuned case is suppressed compared with the resonant 
case, and that the high-purity region is narrowed. From the figure, we can see that for 
coherent photon emission, dephasing should be reduced as much as possible, and that 
for a fixed pure dephasing, k should be taken around n = (2\/2g 2 / e) l l 2 . In addition 
to the main high-purity region, the purity has a small peak (denoted with 'Line 2') 
similarly as the resonant case. From perturbation calculation, we can show that the 
purity has a maximum value of "P max = 3 — 2\/2 ~ 0.17 at 7 P = k/ (2y/2e), 2y/2g 2 /n. 
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(a) On-resonant {u d - u c = 0) (b) Off-resonant (u d - uj c = 8g) 




000 



Fig. 10. The effective purity V(T) and the efficiency VdiT) of photons emitted to the output port 
during < t < T is shown as a function of T. (a) The on-resonant case (w^ = uj c ) for K/g = 2 and 
7p/ff — 0-5, and (b) the detuned case (cj^ — w c = 85) for n/g = \/2/50 and "f P /g = 100. The time T ± / 2 
giving a half efficiency (^(TWa) = 1/2) is also shown by a green line. The unit of time is given as 
T g = K/g, and is about 25ps for g — 25/ieV. 



4-4 Time filtering 

A strategy to improve the purity of emitted photons in to use photons emitted at 
the early stage of decay, because such photons are subject to environment only for a 
short time. The density matrix of photons filtered in the temporal region < t — r < T 
is written as pr{f, r ', t) = p(r, r' ,t)Q(T — t+r)Q(T — t + r'). From pr{r, r' , t), the purity 
of filtered photons are given as 

V{T) = ^^ = i-/' drf dr'\p{r,r',t)\\ (50) 

V d {T) = [Tr(p T )f = (£ T dr p(r,r,t)\ , (51) 

where V d {T) is a square of a probability that a single photon emitted during < t < T. 
V d (T) is square of the probability to obtain a single photon after time filtering. By 
time filtering, the effective purity V(T) is improved at the expense of the square of the 
probability VdiT). 

We first discuss the resonant (u d = uj c ) case. In Fig. 10 (a), we show T-dependence 
of V(T) and VdiT) for n/g = 2 and j p / g = 0.5. In the T — > 00 limit, V(T) approaches 
0.61, which is the purity without time filtering. As T is shortened, the purity V(T) 
increases, whereas the square the probability Vd{T) decreases. If we allow reduction of 
the the square of the probability down to 1/2, one can set filtering time at T = 2r g 
(see a vertical dashed line in Fig. 10 (a)), and obtain the improved purity of 0.85. 
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Improvement of the purity is also possible in the detuned cases. In Fig. 10 (b), we show 
T-dependence of V{T) and V d {T) for n/g = V2/50, %/g = 100, and u d - u c = 8g, 
which is on the upper line of 'Line 2' in Fig. 9 (b). If we allow reduction of the square 
of the probability down to 1/2, the purity is improved from 0.17 to 0.28. 

We can improve the purity further by choosing a shorter filtering time. The efficiency, 
however, decreases exponentially in the limit T — > 0. 

5. Summary 

We analyzed radiative decay of an excited dot in a solid-state cavity QED system 
and investigated the quantum-mechanical properties of the emitted photon, stressing 
the effects of pure dephasing. Our analysis is based on a model in which all the elements 
of the system (including environmental ones) are treated as active quantum-mechanical 
degrees of freedom. We rigorously solved the time evolution of the overall system, and 
derived analytical expressions for the density matrix, pulse shape, spectrum, and purity 
of the emitted photon. These analytical results were visualized under realistic parame- 
ters. The main results are summarized as follows: (i) Changes in the dot decay rate due 
to pure dephasing can be explained in terms of quantum Zeno and anti-Zeno effects. 
The emitted photon pulse length is approximately given by the inverse of the dot decay 
rate, (ii) The emitted photon spectrum agrees with the Glauber formula. The mean 
energy of an emitted photon is not necessarily identical to that of the dot and energy 
conservation is seemingly broken. However, the present analysis revealed that the dot 
can exchange energy with the environment through pure dephasing coupling. Energy 
conservation holds when the environmental energy is included, (iii) The purity of the 
emitted photon is calculated as a measure of indistinguishability. To generate pure sin- 
gle photons, dephasing should be reduced as much as possible and the escape rate of 
cavity photons k should be taken near the optimal value (k; = 2g for resonant case and 
k = (2\/2g 2 /e) 1 / 2 for largely detuned case), (iv) Time filtering improves the purity of 
the emitted photon at the expense of a smaller probability of single photon generation. 

In our approach, we have assumed a white noise for level fluctuation of the QD. 
For detailed comparison with experiments, we have to consider more realistic models 
taking into account phonons 36 ' 37 and background carriers. 25 Extension of the present 
approach to these extended models is an important problem left for future study. 
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Appendix A: Relation to the Master equation and perturbation method 

For analysis on properties of emitted photons, one can employ the Master equation 
approach as an alternative method. The Master equation is derived for the present 
model as 

/ n \ / 'v in —in \ f n \ 

-ig 
W 



d 
dt 



Paa 
Paa 
Paa 
\ Paa ) 



7 
-ig 



ig 

^ + 7 P + iAu 



-ig 




V 









-ig 



2±^ + lp -iAu 
ig 



K 



J 



Paa 
Paa 
Paa 
\ Paa ) 



p(r,r',t) = 



(A-2) 



(AT) 

where p aa = (crV), p aa = (aV), p aa = (o" f a), p aa = (a^a), and Au = u d - u c . Using 
the solution of the master equation, the density matrix p(r, r', t) of emitted photos can 
be written as 

paa(t -r')f3 (r' -r) + p aa (t -r')fi (r' -r) (r<r'<t), 
paa{t - r)(3 (r - r') + p aa {t - r)(3 (r - r') (r' < r < t). 

Here, we have defined (ao(t), (3o(t)) as the solutions of eq. 14 with the initial condition 

(a (0),A,(0)) = (0,l). 

The four eigenvalues of the matrix in eq. (A-l) correspond to /ij's. Therefore, numer- 
ical calculation of eigenvalues of this matrix provides a convenient method for obtaining 
Pi's. It is also useful to treat the matrix in eq. (A-l) in perturbation calculation given 
in Appendix B. Because the matrix in eq. (A-l) is a normal matrix, it is diagonalized 
by a unitary matrix. Therefore, perturbation method used in quantum mechanics can 
be applied to obtain small shift of /ij's with respect to small parameters. 

Appendix B: Expressions of purity in limiting cases 

Simple expression of V can be obtained in some limiting cases. We first discuss the 
largely detuned case {\uJd — ou c \ ^> g). The purity can then be expressed by a sum of 
two contributions as V = Vdot +^cav, where Pdot and V CSN is a purity of photons whose 
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frequency are ojd and u> c , respectively. 35 The former contribution Vdot becomes large 
only below 'Line 1' in Fig. 9 (b), whereas V cstv has a maximum on 'Line 2' in the figure. 
For 7 P , k <C g, perturbation calculation gives 

V dot = — — - x . B-l 

27 P + k 27 P + ne 

The first factor is a weight of a spectrum peak at k = Ud (see eq. (45)), and the second 
one is due to spectrum broadening by pure dephasing. For large detuning (e <C 1), 
Vdot takes 1/2 around 7 P = ek/2, which corresponds to the left-side line of 'Line 1' 
in Fig. 9 (b). For 7 P C g C k, perturbation calculation with respect to g is effective. 
Then, one can consider effective continuum constructed by mixing between the cavity 
mode and the output mode (6-field), and can define a photon emission rate T from a 
dot into the effective continuum. By using Fermi's golden rule, the photo emission rate 
is estimated as T = 2Tcg 2 D(ujd), where D(cj d ) = (K/-n) / ((ujd — oo c ) 2 + (k/2) 2 ) is a density 
of states of the effective continuum. In the presence of pure dephasing, the purity of 
photons emitted from a dot with the rate V is calculated as Vdot = 1/(1 + 27 P /T). Thus, 
we finally obtain 

2q 2 

V dot = — f , (B-2) 

2g 2 + %k 

for k ^> \oJd~ oj c \. One can see that the purity takes 1/2 at 7 P = 2g 2 / k, which corresponds 
to the right-side line of 'Line 1' in Fig. 9 (b). For j p , k <^ g, perturbation calculation of 
V can gives 

2 ^ £K (B-3) 



(/t + 27 p e)(/t + 47 P £:) 

From this expression, it is proved that the purity takes a maximum value 3-2^ 
at 7 P = K/(2y/2e), which corresponds to the lower line of 'Line 2' in Fig. 9 (b). For 
k g 7p; dynamics of the cavity-dot system becomes completely incoherent, and 
can be described by the rate equations 

P**(t) = -T' Paa {t) + T' Paa {t), (B-4) 

Paaif) = +T' Paa {t)-T' Paa {t)-K Paa {t), (B-5) 

where p aa {t) = (a\t)a{t)), Paa (t) = (a t (t)a(t)), and V = 2g 2 /^ p is a transition rate 
calculated from Fermi's golden rule. By solving these rate equations, we obtain 

Vcav = ( K + 2T')( K + ry (B ' 6) 
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From this expression, it is proved that the purity takes a maximum value 3 — 2\[2 at 

7 P = 2\ / 2g 2 /K, which corresponds to the upper line of 'Line 2' in Fig. 9 (b). 

Next, we consider the resonant case (uu c = oud). Some of the features are common 

as the largely detuned case; the purity is reduced for k, > 2g 2 / / ~f p (the right-side line of 

'Line 1' in Fig. 9 (a)), and takes a maximun value 3 — V2 at 7 P = 2^/2 g 2 / k ('Line 2' in 

Fig. 9 (a)). For j p , k <C g, the purity V is approximately calculated as 

= k(2k + 7p ) 

2(k + 1p )(k + 2 1p Y y ' 

From this expression, it is proved that V > 1/2 is realized in the condition k > (1 + 

a/3)7 p ~ 2.737 P , which corresponds to the left-side line of 'Line 1' in Fig. 9 (a). 

Appendix C: Relation between coincidence probability and purity 

Here, we derive the relation P co = (1 — V)/2 between the coincidence probability 
P co and the purity V . We consider the following situation (see Fig. 4). Two solid-state 
emitters (Si and S2) simultaneously and deterministically emit single photons (assuming 
7 = for simplicity) into two input ports (Pi and P 2 ) of a beamsplitter. These two 
photons are mixed by the beamsplitter and are output to ports P 3 and P 4 . We denote 
the photon field operators for port Pj by b(jy (j = 1, • • • ,4), and the pure-dephasing 
reservoir operators for emitter Sj by d^y (j = 1,2). Assuming t — > 00 and 7 = (and 

thus 5 m = 0) in eq. (13), the state vector of Pi photon is given by 

00 „ 

|^i) = J2 drd m x lm {t^x)b\ l)r \ X ), (C-1) 

where \x) = 3^ Xl • • 'd^^ Xm \0). Thus, the emitted photon is entangled with the envi- 
ronment of its source. The input state vector including both Pi and P 2 photons is then 
given by 

00 „ 

1^12)= J2 / drdr'd m xd n x' lm (t,r,xh n (t,r\x')^ iy ^ 2y \x;x'), (C-2) 

m,n=0 

where \x; x') = dj 1)xi • • • $ {1)Xm ^ 2)x , ■ ■ • dj 2)x j0). 

The beamsplitter mixes two photons as b^ r — > [b^ 3 y + b^ r ]/y/2 and b^ 2 y — > [b^ 3 y — 
b^ r }/\/2, but it obviously does not affect the environmental degrees of freedom. The 
output state vector is then given by \ip ou t) = | ^33) + IV'44) + IV'34)) where 

l« = E / drdr'd-x<rx' lm ^ r ' y ' ^ X \ y \ y \x- x'), (C-3) 

m,n=0 
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|« = " E / drdr>d™xd n x> ^ r ' ^ ^ 1»; *'), (C-4) 

m,n=0 

|^34) =2^ d Xd X 2 b (3)r b (A)r'\ X ' X ) 

m,ra=0 

(C-5) 

It is readily confirmed that (^33^33) + (^44^44) + (^34^34) = 1 and (^33^33) = 
(^44 1 ^44 ) • The coincidence probability P co (i.e., the probability of finding single photons 
at both ports 3 and 4), is given by P co = (^34 1 -034 ) • Since the density matrix element 
is given by p(r,r',t) = Ylm=o I ^"^m^'^j^m^,^^), Pco is recast in the following 
form: 

Pu = \ ~ \ j drdr'p(r, r', t)p(r', r, t) = (C-6) 

When pure dephasing is present, the purity of the emitted photon will be less than 
unity and the coincidence probability will be nonzero. 

Appendix D: Proof of eqs. (48) and (49) 

Here, we analytically evaluate E p — j dk k(b' k (t)bk(t))i and E e — j dk k(d\{t)dk(t))i 
in the t — > 00 limit, where (• ■ ■) % = {ipi\ ■ ■ ■ Switching to real-space representations, 
we have 

Ep = \j dr ( {dr ^ r -%( d $))., (D-l) 
E e ='-Jdr {{d r dl)d r - dj(drdj)) . • (D-2) 

From eq. (8), we obtain {{d r b\)b r )i = —K0(t — t)9{t) (( jW)a)j, where r — t — r. Using 
similar equations, we have 

/oo 
dr((^)a-a^(£a)) t , (D-3) 

POO 

E e = -H P / dr - ^(^V)). . (D-4) 

Jo 

Using eqs. (11) and (12) and their conjugates, we have 

r°° ok f°° 

/ dT ^ a) + » dr((a^a) + (^a)), (D-5) 
Jo 2 J 

POO 

E e = gi P dr((a^a) + (^a)). (D-6) 
Jo 

Thus, we need to evaluate I\ = dr(a jf a), I 2 = / °° dr(a j! a) , and I 3 = dr(a j! a) . 
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The equations of motion for (crV), (a^a), and (cr^a) are given by 

l(a^a) = -7(^(7) - ^((^a) - c.c), (D-7) 

— (a^a) = —K(a)a) + ig((a^a) — c.c), (D-8) 

C//6 

j^a) = t(u* d - u c )(a^a) + tg((a^a) - (Ja)). (D-9) 

Integrating these equations with respect to r, we obtain 1 = 7^ +ig{I% — = kI 2 — 
ig(I 3 — I^), = i(oo d — uj c )Is+ig(I 2 — h). When 7 = 0, these equations are solved to yield 
I 2 = 1/k and I 3 = (i/g) x (u*-u d ) / (u c +uj d -u*-uj* d ). Since E p = Ku c I 2 +(gK/2)(I 3 +I£) 
and E e = g^ p {h + ^3), w e obtain eqs. (48) and (49). 
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